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In the arena of the discrete phase space and continuous time, the theory 
of S-marix is formulated. In the special case of Quantum-Electrodynamics 
(QED), the Feynman rules are precisely developed. These rules in the four- 
momentum turn out to be identical to the usual QED, except for the vertex 
function. The new vertex function is given by an infinite series which can only 
be treated in an asymptotic approximation at the present time. Preliminary 
approximations prove that the second order self-energies of a fermion and 
a photon in the discrete model have convergent improper integrals. In the 
final section, a sharper asymptotic analysis is employed. It is proved that in 
case the number of external photon or fermion lines is at least one, then the 
S-matrix elements converge in all orders. Moreover, there are no infra-red 
divergences in this formulation. 

1. Introduction 

It is well known that the Quantum-Electrodynamic (QED) has been 
verified to an extraordinary degree of precision by many experiments. How- 
ever, it can never explain in its renormalized version the theoretical values 
of charge, mass etc. Moreover, it cannot produce in the classical limit the 
exact, non-perturbative solution of the coupled Maxwell-Dirac equations ^. 
Therefore, it is worthwhile to investigate the S-matrix, QED etc. in the dis- 
crete phase space and continuous time which involve non-singular Green's 
functions ^ , in order to obtain convergent results. 

In the section-II, we compare and contrast the continuous and dis- 
crete methods involving perturbation series and Green's functions for ex- 
tremely simple mathematical problems. We conclude that in linear prob- 
lems, both methods incur similar divergence difficulties. However, in non- 
linear problems, the continuous methods encounter divergences whereas the 
discrete methods avoid divergences in the lower order terms. Furthermore, 
the classical potential energy for two particles, as well as the Green's func- 
tions are non-singular in the discrete model. Discussions of this section offer 



suggestive insights into the comphcated problems of QED in subsequent sec- 
tions. 

Wc assume that the free electromagnetic fields and Dirac fields obey 
the difference- differential equations of the paper-II. In the section-Ill, the 
interaction picture is introduced and the relativistic S-matrix is derived in 
the setting of the discrete phase space and continuous time. 

In the next section, we introduce the trilinear interaction energy 
density of QED in the discrete model. We carefully derive Feynman rules 
in the four-momentum space following the difference-differential equations. 
The four-momentum variables do not have cut off limits in spite the presence 
of a characteristic length in the theory. The only difference with the usual 
continuous model is the appearance of the vertex function 5* instead of the 
usual delta functions S. The distribution function 5* is both the joy and 
the pain of this research project. It is a joy since 5* is different from 6 (see 
the Appendix) and there is a slim chance of obtaining the divergence-free 
S-matrix. On the other hand, it is a pain because of the exceedingly difficult 
mathematics involved. We have obtained a simple asymptotic approximation 
d"^ for 5* and apply it to the computations of the next section. 

In the section-V, we apply the Feynman rules for the discrete phase- 
space to evaluate the second order fermionic as well as the photon self energy. 
In both cases, the improper integrals converge. Since the second order terms 
are the most dominant terms in the corresponding S-matrix series, these 
convergences are welcome news. In section- VI, the asymptotic analysis of 
the distribution function 5* and the infinite series representing an S-matrix 
element is sharpened considerably. We investigated an S-matrix element 
corresponding to the physical process involving Eb external boson lines and 
Ep external fermion lines. We prove that in the framework of the discrete 
phase space, the S-matrix elements converge in all orders provided Eb + 



In the present formalism, the basic physical frame-work remains un- 
changed. Only the representations of the quantum mechanics and the quan- 
tum field theory are altered. Moreover, there are no infra-red divergences. 

2. Comparison of continuous and discrete methods in simple prob- 



Let us consider a continuous wave function exp{ikx) in a one dimen- 
sional space. It behaves as the following: 



{3/2)Ef > 1. 



lems 



lim exp{ikx) is not defined; 



x—>oo 



k—*oo 



lim exp{ikx) is not defined; 




(1) 
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Now consider a momentum wave function in the discrete model (see 
equation II-A.I.6): 



{trexp{-e/2)H^{k) 



7ri/42"/2vn! 
n G {0, 1,2,3,...} =: N. 



(2) 



Here, Hn{k) is a Hermite polynomial. In contrast to ([2]), the wave function 
^n{k) satisfy: 



lim ^n{k) = hm C.n{k) = 



k—>oo 



lim / \^^{k)fdk = l. 



(3) 



Now let us try to solve an extremely simple, first order, linear, non- 
homogeneous ordinary differential equation 



ip'{x) = (f){x) 



X e R, 



(4) 



by the method of Green's function. Here, is a prescribed differentiable 
complex-valued function over R. 

The Green's function for (jl]) is given by 



/OD 
k~^exp [ik{x - x)] 
-oo 



dk 



(5) 



1 (x — x) 



2 \(x-x) 



for X X, 



for x = x; 



-i-G{x — x) = 5{x — x) 



lim^^^ I ^G{x - x) I ^ oo. (6) 
The general solution of (jl]), with help of ([6]), is furnished by 



i){x) = a+(C.P.V.) / G{x - x)(p{x) dx 



a -\ — lim 

2 L^oo 



(j){x)dx— I (j){x) dx 



(7) 



where a is an arbitrary complex constant of integration. 
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Now we shall solve the corresponding complex, linear, nonhomoge- 
neous difference equation 



n] 



1 

n G N. 



\/n + lipin + 1) — \/nilj{n — 1) = (t){n)^ 



(8) 



The Green's function for this problem is furnished by 



n) := (C.P.V.) / {kr^Uk)Uk) dk, 

A* G{n,h) = Snh, 
A* G{n,n)\n=n = 1- 

The general solution of ([8]) is given by 

ij{n) = aUO) - ^ (C.P.V.) / k-'Uk)Mk)dk 

n=0 L 



(9) 



}, (10) 



where a is an arbitrary constant. 

Let us apply the solutions as well as ( fTOl) to the momentum 
eigenfunction problem in quantum mechanics, namely 



Pip = eip. 



(11) 



Here, we assume that the parameter e 7^ is sufficiently small. In the 
continuous (Schroedinger) representation of quantum mechanics (ITTi) yields 
the differential equation 



ip'{x) = ieip{x), X E R, ip{x) 7^ 0. 
The general solution of the above equation is obviously given by 



ip(x) = aexp(iex) = a 



1 + iex + 



[lex] 



+ 0(e=^). 



[12) 



(13) 



Here, a is an arbitrary, non-zero, complex constant. 

Now we try to solve (fT2l) by the perturbative expansion and the 
method of Green's function in ([71). Substituting an expansion 



j=0 



(14) 



into ( jT2l) . we derive an infinite string of differential equations: 

V'o(x) =0,^i(x) = ^Po{x),'iIj2{x) =^pi{x),..,^pj{x) = V'j_i(a;), ... . (15) 
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Solving the first equation we obtain 



(16) 



where a 7^ is otherwise arbitrary. In solutions of other equations, we shall 
ignore arbitrary constants. Solutions of the next two equations in (fT5|) by 
the method of (iTl) are listed below. 



V^i(x) 



a 
2 

ax, 
a 

2 

a 



lim 



lim 



dx 



dx 



oc doc 



2'^ 



lim (L^ 



.L—^oc 



(17) 



In the second order solution ip2{x) we encounter a divergent term! Ignoring 
this term ( "renormalizing" ) , by (fT6|) and (ITTj) we can recover the first three 
terms of the expansion in the RHS of (fT3|l . 

Now, we shall try to solve the same problem fllip with the discrete 
representation of quantum mechanics. The corresponding difference equation 
and the exact general solution are provided by 



A*'0(n) = ieil>{n) 
ijj^n) = a^n{e) = a 
a ^ 0. 

Using a perturbative expansion 



ao)+eC(o) + -e::(o) 



+ 0{e'), 



j=0 



(18) 



(19) 



the difference equation in (fTSll yields the following infinite string of difference 
equations: 

A*Mn) = 0, A*Mn) = Mn), A*^j{n) = ^j-i{n), ... (20) 
The first of the equations (!20|) is solved by (see Appendix-I of paper-II) 



^o(2n + l) 



a^ri(O), 
a6n(0), 
aWilO) = 0. 



(21) 
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The second and third equations of (!20|) can be solved using (fTOjl and (12T|) . 
Ignoring arbitrary constants, these solutions are: 



^i(2n+ 1) 



-la 



/oo 
k'^in{k)b{k) dk 
-oo 



= 0, 



-ta{ lim A;-^ [e2n+i(A;) - WilO)]} 

fc-^oo 
-^<2n+l(0), 

-a{(C.P.V) 



-a 



eL(0) - hm [A;-26n(A;)] 



E^2n+l(0)6n+l(fc) 



.n=0 



7/'2(2n + l) = 0. (22) 
(In deriving equations (122|) . we have used 

oo «QQ 

E ^2n+i(0) = -^\k), / /(fc)5'(A;) dA; = -f\k) (23) 

etc.) The RHS of the equation for %l)2{2n) has a divergent constant. Ignoring 
it and using equations f|T9|) . fl2Tl) and fl22l) we can recover the three terms of 
the RHS of the series in equation f[T5]) . So, we discover that the solution of 
the linear quantum mechanical problem in flTT]) with perturbative expansion 
and Green's functions produces exactly similar divergence difficulties in the 
continuous or the discrete representation. 

Now let us investigate a simple non-linear toy model by the two 
methods. Consider a first order non-linear ordinary differential equation: 



The exact general solution is given by 

iplx) = a[l — eaxY^ 



(24) 



(25) 



where a is an arbitrary complex constant. In the case where a is a non-zero 
real constant, the solution in fl2^ has a singularity at x = (ea)~^. In case 
\ax\ < 1, we can elicit from (125!) a series expansion 



a ^^(eax) 

j=0 



(26) 



A perturbative expansion 



^(^) = X^(e)^V'j(a;) 

j=0 



(27) 
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leads to the following string of differential equations: 

= 0,V^;(x) = [ij^{x)]\i,'^{x) = 2^oWV^i(x),... . (28) 

Using the solution ([7]) involving the Green's function and ignoring arbitrary 
constants after the first solution, we obtain 



a lim 

L^oo 



xdx — xdx 

L Jx 



a^x^ - hm L^. 



(29) 



Ignoring the divergent term in the last equation, using equations (1271) and 
( |29|) . we recover the first three terms in the expansion of (126|) . 

The corresponding non-linear difference equation is furnished by 



With a perturbative expansion 

oo 
j=0 



(30) 



(31) 



the equation ( ]30l) implies that 

A*Mn) = 0, A#^i(n) = [Mn)]^A*Mn) = 2^o(n)^i(n), ... :: . (32) 

The solutions of these equations by the method of Green's function in (fTOj) 
lead to the following expressions: 



^i(2n) 



ae«(o), 



^ r poo 

$^{[6n(0)]2 (Cpy^) / k-%n{k)'Uik)dk 
n=0 I J-oc 



(33) 



^i(2n + l) 



/oo 
{A;-i6n+i(A;) 



-ia^ (C.P.V. 



oo 
oo 



.n=0 



(27rfc)-^6n(A;)[4^5(A;) + 



(34) 



(Here we have used A. 18.) In the last integral there are neither "ultraviolet" 
nor "infrared" divergences. 
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The third equation in (|32]) yields the solution 



^ ^ poo poo 

Mn) = -2a=^EE{[^2.(0)][e„(0)] (C.P.V.) / / {kk)-' 

^2h{k)£,2n{k)in{k)i2h{k) dk dk } 
= 0. (35) 

In the RHS of equation fl5^ . the (C.P.V.) integral converges to zero since the 
integrand is odd with respect to the variable k . Therefore we conclude that 
for a non-linear problem, convergence is more likely in the discrete method 
compared to the continuous case. 

No we shall discuss the vertex function in the one dimensional mo- 
mentum space in the case of a trilinear interaction term. In the usual con- 
tinuous case, it is given by 



/OO 
exp [i{p-q-k) x] dx = 5{p- q- k). 
-OO 



(36) 



The above equation indicates a sharp conservation of momenta or wave num- 
bers for three particles. In case functions are defined in the half line x > , 
the corresponding vertex function is furnished by^ 

(27r)-i exp [t{p-q-k) x] dx = ^{5ip-q~k) + '-'-i^^^^^}. (37) 

In the RHS, there is one term indicating the sharp conservation and another is 
representing one soft conservation. In the field theory in an one dimensional 
lattice space, the corresponding vertex function is^ 

OO OO 

{27t)-^ exp[i{p - q - k)n] = ^ 6{p - q - k + 27cj). (38) 



J=-CO 



The above indicates a denumerably infinite sharp conservation of wave num- 
bers. In a single discrete phase plane, the corresponding vertex function is 
specified by (see equations A.l): 

OO 

6*{p,-q,-k) := ^a(p)^„(g)en(fc) ^d*{p,-q,-k). 

n=0 

In the special case of p = 0, it is exphcitly furnished by (see A. 18): 
d#(0,-g,-A;) = -{^\^Tk\5{q + k) + ^\^^\5{q-k) 

TT 

+ -A\q + k\-'/' + \q-kl^'']}. (39) 
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The RHS of the above equation demonstrates two sharp as well as two soft 
conservations. Moreover, spontaneous reflections of momenta (on lattice 
points) are allowed. (These reflections are analogous to Bragg reflections 
on the lattice planes of a crystal.) 

Now we investigate the potential equation in our discrete model. It 
is given by: 



5'^^AfAfy(n) = 0, (n) := {n\n^,n^) G 



(40) 



The corresponding Green's function (which is the potential energy between 
two unit charges) satisfy: 



G{n,h) = 
lim \G{n, n)| < oo 



fk-k) 



a=l 



(41) 



The convergences of the integrals above are due to the facts that Hna{ka) is 
a polynomial and there is a decaying weight factor exp[— (k ■ k)] in products 
of in'^ika). 

Now we consider the Green's functions (equation (II-A.II.l.B)) for 
the difference-differential Klein-Gordon equation. The coincidence limits are 
provided by: 



A(„)(n,t;n,t;/i) = (27r) 



TT 



'3/2 



{exp[-fc^ -kl- kl] 



n 



\ 2"'(nJ')! 



{v'"^k^kp + fi^) ^ dk^ 



(a) 



}rf^k. 



(42) 



Because of the weight factor exp[— A;^ — k^ — the integral on the RHS 
of (H2|) converges. All of the Green's functions for free fields in our discrete 
model are non-singular. 

The discussions of this present section will provide valuable insights 
into the complicated topics of the following sections. 

3. The interaction picture and the S-matrix 

We shall follow the same notations as in the previous papers I and II. 
The equations of the relativistic quantum fields are expressed exclusively by 
the partial difference- differential equations in this paper. We shall now derive 
the S-matrix for interacting fields in the sequel. In the interaction picture, 
the time-evolution of a Hilbert space (state) vector \^i{t)) (representing a 
many particle system) is governed by the differential equation^ 



idt\^i{t)) =Hi{t)\mi{t)) 



(43) 
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Here, Hi{t) represents the hermitian operator corresponding to the total 
interaction energy at the instant t. We can express Hj{t) in terms of the 
interaction energy density operator 7i/(n, t) by the triple sum: 

oo (3) 

Hj{t):= Yl ^Kn,t), (44) 

n=0 

Hiit) = -Ci{n,t). (45) 

Here, Cj{n,t) stands for the interaction Lagrangian density and it is a rel- 
ativistic invariant operator. (The equation ( l45l) holds for most of the useful 
cases.) A necessary micro-causality requirement is 

[niin,t),niin,t)]=0 forn^n. (46) 

Usually, the interaction energy density operator 7i/(n, t) has an overall mul- 
tiplier which is small. It is customary to solve ( l43l) by perturbative series in- 
volving this small parameter. Moreover, the S-matrix is the operator which 
takes a prescribed initial state \i) into another prescribed final state |/) con- 
sistent with the evolution equation fHSl) . The perturbative series expansion 
for the S-matrix is furnished by^ : 

oo oo oo (3) oo (3) 

S = / + 1]5,=/ + I][(-0Vj!]E / / ^^2... 

j=l j=l ni=0 n2=0 

oo (3) 

J2 / dt,{T[nj{n,,t,)nj{n,,t,)...n,{n„t,)]}. (47) 



Here, T stands for Wick's time ordering operator. There is another operator 
N called normal ordering^. It arranges creation operators to the left of 
annihilation operators. Furthermore, there is still another operation called 
contraction between two operators^ and it is defined by: 

A{n, t)B{n, t) := T[A{n, t)B{h, t)] - N[A{n, t)B{h, t)]. (48) 

We can extract from the commutation and anti-commutation relations (II- 
3.10B), (II-4.9B), (II-5.11Bi-vii) and the linear relationships (II-A.II.4B), (II- 
A.II.Sc), the following examples of contractions: 

A^{n,t)A^{h,i) = -ir]^^DF{n,t;h,i)I, (49) 

^(n,t)^(n,t) = ^(n,t)§(n,t) = 0, 

^(n,t)^f(n,t) = -'^{h,i)'i!{n,t) = iSF{n,t;n,i;m) I, 

^(n,t)A^(n,t) = ^{n,t)A^{n,i)=0. 

The Green's functions Af{..) , Df{..) , and Sf{..) are all defined in Appendix- 
II of paper-II. (We should mention that these Green's functions are analogues 
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of the causal Green's functions of Stuckelberg in the +2 signature and differ 
from the corresponding Feynman- Dyson propagators by a factor (— ^2).) 

Now, we shall state the Wick's theorem'^ on the decomposition of a 
chronological product in the arena of the discrete phase space and continuous 
time. It can be succinctly stated as: 

T[A(ni, ti)5(n2, t2)C(n3, h) JK, t,)] 

= N[A{m,ti)B{n2,t2)Cins,h) J(n,-,t,)] 

+N{[A{ni,ti)B{n2,t2)C{ns,t3) J{nj,tj)] + all other single contractions} 

+.... 

+N{[A{ni,ti)B{n2,t2)C{ns,t3) J{nj,tj)] + all other double contractions} 

+.... (50) 
+A^{terms with maximal number of contractions}. 

(This theorem is provable by induction.) 

4. Feynman rules of Q.E.D. in discrete phase-space and continu- 
ous time 



We choose the interaction density as 

ni{n,t) := -teN[^{n,t)Y'^{n,t)A^{n,t)]. (51) 

Here, the parameter |e| = y^Air/lSY is a small positive number. The in- 
teraction term in (15T1) is derived from the principle of minimal electromag- 
netic interaction. According to this principle, the difference and differential 
operators A*, dt in the free Lagrangian density (II-5.3) are replaced by 
A* — ieAa{n,t) and dt — ieA^^Uji) respectively. Moreover, from the discus- 
sions in section-IV of paper-I, it is amply clear that the interaction energy 
density in flSTj) is a relativistic invariant. Furthermore, the micro-causality 
condition (3.3) is satisfied by interaction flSTl) due to equations (II-3.10B), 
(II-4.9B), and (II-5.11B i-vii). 

The equations (3.4) and (IHTj) yield for the S-matrix: 

oo oo oo (3) oo (3) „ 

5 = / + 5^5,=/ + 5][(-e)Vj!]$^ / dt,...J2 / 

j=l j=l ni=0 nj=0 

T{N[^{n,,t,)Y''^{ni, t^)A^, (m, t^)]... 

N[^inj, t,)Y^^in„ tj)A,^ (n„ t,)]}. (52) 

The R.H.S. of fl52|) is a relativistic invariant operator. Let us consider the 
second order term 5*2 in fl52l) for the sake of simplicity. By the Wick's de- 
composition (I5U]) and equations (H^-v). 5*2 can be reduced to the following 
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sum^ : 



52 = (eV2)$^ \ dt^Y. dt2[J2^iA){ni,h;n2,t2)] (53) 

ni=0 "'^ 02=0 "'f^ A=l 

8 

= : YS2(A)- 



A=l 



Two of the operators r(^)(..) are explicitly furnished below as: 

r(5)(ni, ti; 112, h) = -A^[7^^(ni, ti)S'F(n2, ni, ti; m)T]^^DF{n2, ^2; ni, ti)^{n2, ^2 
r(6)(ni,ti; n2,t2) = r(5)(n2, t2; ni, ti). 

The second order terms 82(5) = 5*2(6) contribute towards the self-energy of 
an electron. 

Let us now work out the matrix entry (/|r(5)(..)|i) for the initial 
state of one electron \i) := aJ(p)|\l/o) and the final state of one electron 
I/) :=at(p)|vl/o). 

Using equations (II-5.7Bii), (II-5.8), (II-5.10), (II-A.II.5B), and (II- 
A.II.6B), the operator r(5)(..) in ( 15^ yields 

(/|r(5)(..)K) = m[E(p)E(p)]-V2a^(p)^M{iim / / d'kdS 



.b=l 



(27r) exp[t{E + ki + ^4)^2]^?, 



(55) 



a=l 

-- d^qdq'^ = —d^qdqi. 



(27r)~^ exp[-i{E + h + g4)ti]}u^(p) 



We shall now define a new distribution function (see the Appendix) 
by the triple sum 

00 (3) 3 

(p, q, k) := {[[ U {pj)U (%)] 

n=0 j=l 

00 00 00 

ni=0 rfi=0 n3=0 
[^2 (P2 ) ^2 ( 92 ) ( ^2 ) ] [^n3 (P3 ) ^3 ( 93 ) ^3 ( ^3 ) ] • 



(56) 



Therefore, by the equations (l53l) . (l56l) . and ( ]56l) . we obtain the second order 
contributions for a fermionic self-energy as 



(/|52K) = (/|52(5) + 52(6)K) = (2!)(/|52(5)K) 
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5iE - - A;^)7^(gag" + - ity\iYqp - ml)r]^,{k^k(, - te)-'Y (57) 

53#(p,-q, -k)%^ + A;^-i5;)}u,(p). 

A graphic way to represent the right hand side of the equation (1571) 
is by Fig.l. 

(P.E) 




(P.E) 



Figure 1: The second order fermion self-energy graph in the momentum space 

This expression and other matrix elements give rise to table -1 of 
Feynman Rules. 

The Feynman rules in the four-momentum space (instead of the rules 
in the discrete phase space and continuous time) will be used in the sequel 
for the sake of simplicity. Integration of all the internal (or virtual) bosons 
and fermions (or anti-fermions) must be performed over R^. (The fourth 
component of the momentum is integrated over the real axis because of the 
addition of —ie in the propagators.) The limit e ^ 0+ is taken after all the 
internal integrations are performed. 

However, the actual construction of a matrix element < f\Sj\i > 
in (l52l) from Table-I is somewhat incomplete until we determine the correct 
numerical factor which multiplies the ordered product of operators with the 
same j and the same physical process. It can be deduced that the multi- 
plicative factor must be 

Kr.= {-iysgn{a){-q\e\y. (58) 

Here, / is the number of closed loops, a is the permutation of the final 
fermions, and q = —1 for electrons, ±1/3 for quarks etc. The integer j 
is the number of vertices. In case, (5*(p, — k, — q) in Table-I is replaced 
by (27r)^(5'^(p — k — q) , we obtain exactly the same Feynman rules which 
emerge out of the usual theory. However, we shall prove in the Appendix 
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Factor In S-Ma1rix 
{dhscrete phase-space 
& contnuoue time} 


Factor in5;JJ^tt 


1 ucHTipuriQni 


Photon PripAg^m. 










I 
















tfiooirtrit om^lsotrig 

external photon lines: 




yy 


Incoming or outgoing 






m 


Inoumfng oroufgoing 
axtemflianti-fermion lines: 

^ \ . .' ' . 1 
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Table- 1: The second order photon self-energy graph in the momontum space. 



that Sf{p, — k, — q) ^ (27r)^5^(p — k — q) . Feynman rules in Table-I are all 

manifestly relativistic except possibly the vertex term. 

According to the Feynman rules in Table-1, energy is precisely con- 
served at each vertex. Therefore, we can introduce a physically meaningful 
matrix M* by: 

oo 

{f\M#\i) : ^{flJ^Mfli), 

{f\Sj\i) = :ii27r)5{E^f)-E(,)){f\M*\i), 
{f\S-I\i) = i(2n)5(E^f)-E^^){f\M*\i), (59) 

where E(^i) and E(^f) are the initial and final energies respectively. 

The transition probability from the initial state \i) into the final 
state I/) per unit time is provided by: 

^(/)(,) := {27rf 6{E^f) - E^^)\{f\M*\t)\\ (60) 
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5. The second order self-energies of electron and photon 

For the analysis of the second order self-energy of an electron we start 
with the equation (j57l) . Moreover, utilizing (!59l) we can write 

{f\S2\t) = z(27r)5(i - E){f\M*\t) (61) 

# 

= : -2e'm[E{p)E{p)r'/'Mp)J2^P,P)Mp), 

(2) 

-2V(p,p) : = lim / / d^qd^k6*{P, -q, -k)6{E - - k^) 

x53*(p, -q, -k)5(-E + + k^)l^{iYqp - m/) 

X (g-g, + - ie)-\kPkp - ie)~^^^. (62) 

Using the properties of the Dirac matrices and the asymptotic ap- 
proximation (i*(p, q, k) in equation (A. 18), we obtain from (l62l) that 



V(p,p) ^ hm / / Ac^'fc4(p,-q,-k)5(i-g^-fc^)4(p,-q,-k) 

^ ^-^0+ Jr* yR4 

5(-E + g^ + fc^)(z7''gp + 2m/)(g"g„ + - ie)~^ 
x{k^kp-iey^. (63) 

The last two factors in the above integrand can be combined by the parameter 
X into 

{q''qa,+m^-it)-^{k'^kp-it)-^ = [ dx [x(g"g« + m^) + (1 - x)k'^kp - ie] . 

Jo 

(64) 

Moreover, we make the simplifying assumptions 

p = p = 0. (65) 
By dSl-ES]) and (A. 15) we obtain 
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V* (0, m; 0, m) = lim [ [ [ d^qd^kdx\df{0,-q,-k)\^ 
[6{m -q^- k^)f {iYqp + 2m/) [a;(g"g„ + m^) + (1 - x)k'^k,3 - le] 
= (4)-i(47r)-6 lim /" [ [' d^kd^dx [6im - q^ - k^)?Y\ (A^{\ka + qa\ 

^^°+ jRi jRi Jo 

[5{ka + qa)f + \ka " qa\ - qa)f + 2\kl - ql\^/H{k^ + q^k^ - qa)} 

+2{Af{6{ka + qa) + ' qa) + \{ka + qa) / {ka ' qa)\'^^S{ka + qa) 

+ \{ka-qa)/{ka + qa)\'^'S{ka-qa)} 

+ 4 [\ka + qa\-' + \ka - ^al"' + 2\kl - ql\-'/']) [i^qp + 2m/) 

+ m^) + (1 - x)k^kp - ie] . (66) 

To analyze the above integral we employ some mathematical devices, 
(i) We put ^/\k\Kk) = \k\ [5{k)f = 0, ^{{k + q)/ik - q)\6ik + q) = for 
k 7^ —q. (ii) Whenever there is a factor in the integrand which is odd with 
respect to some integration variable, we drop that term. With such tricks, 
the RHS of (l66l) can be reduced to the form 

4 V* (0, m; 0, m) = {A)\Att)-'^6{0) [A^A + A\B^ + B2 + B3) + A''{Ci + C2 + C3) + D] 

(2) 



A := hm [ d^k [ dx [-i'j\m - k^) + 2m/] [k^kb - {k 



4 \2, 22 •l"^ 

mx) + m X — ie\ 



(67) 

V- 

(68) 



5i : 



lim ! dqi [ d^k [ dx [-i-f\m - k^) + 2ml] [k^h-{k 



4 \2 

mx) - 



2 2 
m X 



+ a;(g? - A;?) - ^e\ [|A;i + g^l + \k, - q,\ + 2\kl~ ^2] 



: 



lim / dqidq2 / c^'^A; / dx [— 'i7^(m — A;"^) + 2m/] [A;^A;;, — {k 



mx)^+ 



2 2 
m X 



+ x{ql + g^' - A:? - fc^') - ^ [1^2 + q2\-' + 1^2 - g2r' + 21^2^ " ^ll'"^] 



D : 



lim ! d\ [ d^k [ dx [-i^im - k^) + 2m/] [A;^A;b - (A;^ 



2^2 



mx) + m X 



+x{q^qc — k'^kc) — ie 



-2 



n (l^a + qa\-' + |A:. - g.r^ + 2\kl- ql\-'/') 



a=l 



(69) 
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In the above expression, a mathematically unacceptable symbol, S{0) 
occurs. We can explain it by physical arguments of energy conservation 
6{E -E) = 5(0) in (EH]) and ([59]). 

We note that A, Bi^Ci, D are all "logarithmically divergent" inte- 
grals. Therefore, symbolically we can write: 



'(2) 



m) = m 



(5(0) [log. div. terms] 



(70) 



From the 5*2- matrix element of the usual quantum field theory in the space- 
time continuum, we obtain 

V (0,m;0,m) = m[6{0)f[ d^k f dx^kf^kp + m^x^Y^ [21 - ~ x)] 

= m [5(0)]*^ [log. div. terms] . (71) 

Comparing (ITOjl and (!7T|) we conclude that the 5*2 matrix element for the 
fermionic self-energy case in the discrete phase space and continuous time 
does converge! 

We shall now work out the second order photon self-energy (see fig. 2. 

Below.) 




Figure 2: The second order photon self-energy graph in the momontum space. 



This self energy is given by the Feynman rules in Table-I as propor- 



tional to 



n#M'^(k) : = lim{Tr / / d'^pd'^pY 



i'~f''Pp — ml 
p'^pa + m'^ — 

3 

n l'^*(-Pc,Pc, fee 
c=l 



p^Pp + m? — ie 

2 



r [si-p^+p' + u] 



(72) 
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We use the asymptotic approximation (A-I.21) to evaluate n*'*''(k) 
above. It is given by (dropping e) 



n#M^(k) 



-(47r)-^(5(0)Tr{7^ 



) 



A :-- 



+ {2)-'(Ci + C2 + C3^+{8)-'d'\y}, 
/ dfp [i'y'^Pa — ^I] [i'y^ip'^ — u) + ml] [p^Pf^ + m^] 



(73) 



X 


{ [{pi + hf + (P2 + + (P3 + hf - 


- (P4 - 


-vf + m^^ 


+ 


[{pi + kxf + (p2 + + (ps - k^f - 


(P4 - 


vf + mY^ 


+ 


[{Pi + kif + (p2 - k2f + (P3 + A;3)' - 


(P4 - 


vf + m^^ 


+ 


[(Pi + k,f + (p2 - k2f + (P3 - A;3)' - 


(P4 - 


uf + m^^ 


+ 


[{pi - kif + {p2 + A;2)^ + (P3 + k^f - 


(P4 - 


vf + m^]~^ 


+ 


[{Pi - ki f + (P2 + /C2)' + (P3 - ^3)' - 


(P4 - 


z/)2 + m2]-' 


+ 


[iPi - kif + {P2 - k2f + (P3 + k^f - 


(P4 - 


i/)2 + m2]-' 


+ 


[(Pi - A;i)^ + (P2 - A;2)' + (P3 - ^^3)' - 


(P4 - 


^)^ + mT'}, 



R4 ^R 



ml 



\}l'^{p'^ — i/) + m/] 



X { [ [p? + (P2 + k2f + (P3 + ^3)' - (p' - i^f + m^] 

"p? + (P2 + k2f + (P3 - A:3)' - (p^ - uf + m^] 



+ 
+ 
+ 

X 

+ 
+ 
+ 
+ 
+ 
+ 



+ m^] ^ 



pl + (P2 - k2f + (P3 + k^f - {p^ 
p\ + (P2 - k2f + (P3 - ^3)' - (p' - i^f + m"] 
Pl -Pi - ki\~^ + |pi +pi + ki\~^ + |pi -pi + /cil"^ + |pi +pi - A;i|~^ 
Pl -Pi -A;i)(pi+pi + A;i)r^'^^(l + sgn(pi -pi - A;i)sgn(pi +pi + /ci)) 
Pl -Pi -^i)(pi -pi + A;i)r'^''^(l + sgn(pi -pi - A;i)sgn(pi -pi + ki)) 
Pl -Pi -^i)(Pi+Pi - A;i)|~^''^(l + sgn(pi -pi - A;i)sgn(pi +pi - /ci)) 
Pl +Pi + /si)(pi -Pi + A;i)|"^/^ (1 + sgn(pi +pi + A;i)sgn(pi - pi + ki)) 
P1+P1 + ki){pi +pi - A;i)|"^/2 (1 + sgn(pi +pi + A;i)sgn(pi +pi - /ci)) 
|Pi -Pi + A;i)(pi +pi - /ci)r^/^ (1 + sgn(pi - Pl + A;i)sgn(pi +pi - /ci))]} 



(7 := / d^p / dpi (ip2 [i7"Pa — ml] [i'j^{p'^ — z/) + m/] [p^pp + m^] 

X { [Pi + P2 + (P3 + ^3)' -{p'-i^r + mY'+ [Pl + P2 + (P3 - ^3)' - ip' - + m^] 
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2 

X {Jl [\pa + Pa + kA\~^ + \pa-Pa - /^aI""^ + \pa+Pa - ^aI""^ + \pa - Pa- ^aI" 

A=l 

+ \{Pa+Pa + kA){pA -pA- kAjl"^^"^ (1 + sgn{pA +PA + kA)sgn{pA -pA- kA)) 

+ \{Pa+Pa + kA){pA +pA - kA)\'^''^ (1 + sgn(p^ + + fcA)sgn(pA + Pa - ^a)) 

+ Kpa + Pa + kA){pA -pA + ^A)r^^^ (1 + sgn(p^ + p^ + /cA)sgn(pA - Pa + ^a)) 

+ \{PA-pA - kA){pA + Pa - /i;A)r"^^^ (1 - sgn(pA - Pa - A;A)sgn(pA + Pa - /^a)) 

+ Kpa - Pa - ^a)(pa -pA + kA)t^^'^ (1 - sgn(pA - Pa - A;A)sgn(p^ - Pa + ^a)) 

+ Kpa + Pa- kA){pA - Pa + ^A)r^^^ (1 + sgn(p^ + Pa - /cA)sgn(pA - Pa + ^a))] } 



/ (i'^p / (i'^p [«7"Pq: — mJ] [i7'^(p^ — z/) + mJ] 
[p'^P/3 + m^] [Pbpb - {p'^ -^T + fn^] ^ 

3 

{JJ [be - Pc - kcl^'^ + Ipc + Pc + z^cT^ + Ipc - Pc + /^cT^ + \pc + Pc - z^cT^ 

c=l 

KPc - Pc - /i;c)(Pc - Pc - /i;c)r^^^ (1 + sgn(pc - Pc - /^c) sgn(pc + Pc + /Cc)) 
|(Pc - Pc - /Sc)(Pc - Pc + ^c)r^^^ (1 - sgn(pc - Pc - sgn(pc - Pc + ^c)) 
KPc - Pc - /^c)(Pc + Pc - /i;c)r^^^ (1 - sgn(pc - Pc - /^c) sgn(pc + Pc - A^c)) 
KPc + Pc + kc){pc -pc + kc)\''^^^ (1 - sgn(pc + Pc + fee) sgn(pc - Pc + kc)) 
KPc + Pc + ^c)(Pc + Pc - ^c)r^^^ (1 - sgn(pc + Pc + fee) sgn(pc + Pc - /Cc)) 
KPc - Pc + /i;c)(Pc + Pc - /i^c)!""^^^ (1 + sgn(pc - Pc + /ec) sgn(pc + Pc - A^c))] } 

(74) 

The analysis of the integrands in A, Ba, Ca, D imply that 

jj*/^!^ _ ^(-g-) [quadratically divergent terms] . (75) 

However, the expression from the usual theory in continuous space-time in- 
dicates that 

li^^ = [5(0)]^ [quadratically divergent terms] . (76) 

Comparing fl75|) and fl76l) . we conclude that the second order photon self- 
energy integral converges in the discrete phase space and continuous time 
formulation. 

6. Convergence of the S-matrix element {f\Sj\i) 

We start this section with some simple examples. Let us consider the 
improper Riemann integral k dk . This integral does not converge. How- 
ever, the Cauchy principal value CP. V. k dk "converges" to zero. Now 
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consider another improper integral k"^ dk . It diverges (strongly) and its 
Cauchy principal value does not converge either. Now let us investigate the 
infinite series 1 — a; + — (3!)x^ + (4!)x^ — ... . It is a divergent series 
for a; 7^ 0. However, it is Borel-summable to the value e~^[l + xt]~^dt 

oo 

(which converges). Consider another (complex) infinite series Yl exp[i2"a;]. 

n=0 

It is strongly divergent for all x G R and cannot be "summed" in any sense. 

In this section, our primary goal is to investigate the convergence of the 
j-th order S-matrix element (/ l^jl i) , where \i > and \f > are initial and 
final states respectively. This matrix element is furnished by the equation 
( 152|) . We apply Wick's theorem inherent in ( 150|) and consider only terms 
with a specific number of contractions. There exist many terms even in that 
category. We pick a typical term (^f \Sj \ i) in a particular category. Using 
the Table-I, we can express the matrix element as the following infinite series 
of integrals: 



(const) 



oo(3) 

E 

ni=0 



/ dt,...y2 / dt,{< f\N[r^ij{n,,t,)A,A^U' 
Jr. „^r,JR 



'y^'^Spini, ti, 112, t2] m)?7^2^3Dir(n2, ng, t^)SF{jy2, ^2; ng, tg; m)7''^.. 
Spijcvj-u tj-h nj, tj-; m)ip-f^' Af,^ (n^, tj)]}. (77) 

The right hand side of the above equation contains multiple series and 
multiple integrals. There exist three possibilities. It can be convergent, or 
summable in some sense, or strongly divergent. Before embarking upon such 
investigations, let us consider a similar problem for a much simpler series of 
integrals: 



E 

n=0 



fn{k)dk 



A necessary criterion for convergence is that fn{k)dk converges (or con- 
verges as Cauchy principal value) for each n G {0, 1, 2, ...} to If further- 
more, Ylo'i^nSn) converges, then the original series of integrals is said to 
converge. In case the series Yl'^{cnSn) does not converge, we should deter- 
mine whether it is summable or strongly divergent. (In the case of the right 
hand side of ( 1771) . such analysis is extremely cumbersome.) In case the anal- 
ysis of Yl'o'^^n^n) is very difficult, we can try the asymptotic analysis. In this 



CnjT^fn{k)dk 



is summarized 



strategy, criteria for "convergence" of Yll. 
by the following steps: 

i) Prove that C.P.V. fn{k)dk converges to s„ for every n G {0, 1,2, ...}. 

ii) Express the series as X]^=o '^n jT^ fn{k)dk +Y1'^=n+i jT^o fn{k)dk 
a sufficiently large positive integer A^. 

ni) Prove that the series Yl'^=i^N+nfN+n{k) converges to (Jjq{k) . 



for 
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iv) Check that hm (JN{k) = 0. 

v) Finally, prove that C.P.V. J^^aN{k)dk = [Z]^=Ar+i Cn/n(/c)] dk 
converges. 

We shall follow a similar five step strategy to investigate the convergence 
of the matrix element |S'j| in (177|) . For the step i), we have to inves- 
tigate improper integrals representing propagators (See the Appendix-II of 
the paper-II.) These are furnished by: 



D(,)(n,t;n,£) = {27i)-' [ 



.b=l 



[v'^'^kakp) ^ exp [ik^it - i)] dk^ 



(a) 



ci^k, (78) 



%)(n,t;n,t) = (27r)-i / 



.b=l 



{i'j'^Pfi + ml) ^ exp [ip4^{t — £)] dp'^ 



(a) 



d^p. (79) 



With the help of reference we have integrated explicitly the right hand 
side of fl78l) in case of D+(n, t; n,t). The result is the following: 



D+{n,t-h,t) = -% [2(2v^)3] 'exp {-{t-tfl%\ 

[ni/2] [hi 1 2] 



.6=1 



ii=o j3=o 



ni+ni-2(ji+jfi) 



\ _|_ (^_]_^"2+n2-2(j2+j2) 



xF 



-1 



xF 
xD 



Ji!..j3!(ni-2ji)!..(n3-2j3)! 
2-1 (^m + rii - 2(ji + ji)) + l] F [2-1 + n2 - 2(j2 + J2) j + 1 

(ns + ria - 2(j3 + ja)) + l] (r [2-1(^1 + ... ^n^- 2(ji + .. + J3)) + 3 

-[(ni+...+n3-2(ji+..+j3))+2] (^(^ " ^)/v^) ; 



: = n/2 for n even, 
{n — l)/2 for n odd. 

Here, D-p{iz) is the parabolic cylinder function 1°. 



(81) 
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The right hand side of (IHOl) is defined (or non-singular) everywhere. In 
fact, in the coincident point, 

D+ (0, 0; 0, 0) = -i^/^ = -D_ (0, 0; 0, 0) . (82) 

Thus, by equations (A.II.4B) and (A.II.5B) of the paper-II, Dp (n, t; n, t) in 
(ITS]) is non-singular everywhere. Similarly, the integral representing Sp (n, t; n, t; 
in (17^ converges by the equations 

(n, t; n, t; m) = (7^ A* — mJ) A^? (n, t; fi, t; m) ; 
Dp (n, t; n, t) := /\f (n, t; n, t; O) . (83) 

In fact, in the coincident points, 

5+(0,0;0,0;m) = 2"^ [-7^ + im3^(3/2, 2; m^)/] , 

5_ (0,0;0,0;m) = 2"^ [-7^ - ^771^^(3/2, 2; m^)/] . (84) 

Here, ^'(a,/3; z) is a degenerate hypergeometric function 

In the integrations for (n, t; n, t) in flHOj) and S-t (0, 0; 0, 0; m) in 

(6.8) neither ultraviolet nor infrared divergences are encountered. 

Since wave functions A^{n,t) and ip{w,t) are defined everywhere, 

we can conclude that (initial) finite part of the infinite series in ( 1771) . namely 

Ni(3) ^ N,(3) 

[cons'^ 



t)J2f ^^1- E / dt,{<f\..\t>} (85) 

ni=0 "'^ n,=0"'^ 



is always non-singular. We have to investigate the convergence of the re- 
maining tail-end: 

00(3) 00(3) 
(const) / dti... / rftj{< f\..\i >}; (86) 

ni=Ni+l"^^ nj=Nj+l'^^ 

where 1 := (1, 1, 1). For that purpose we take recourse to the analogue of 
step iii) in the strategy. We have to prove that 



00 „ 00(3) r r r 

y I dti... 2, / dtj / d'^ki... / (i'^pi..! Fourier Transform of < f\. 



00(3) 

E 

ni=Ni+l''"- nj=Nj+l 

(87) 

converges. In this endeavour, we use Fourier transforms of all fields (as in 
(4.5)). Thus we encounter at each vertex the modified distribution function: 



53#v(p,q,k) := J2 J2 Yl 

ni=Ni+l n2=N2+l n3=N3+l 
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.b=l 



(Compare the above with the distribution function (5*(p,q, k) in the Ap- 
pendix.) Let us first try a simpler distribution function: 

oo 

4(p,q,k):= ^ Up)Uq)Uk), {p,q,k)eR'. (89) 

n=N+l 

The exact sum in (!89|) is intractable presently. Therefore, we try the asymp- 
totic approximation for a sufficiently large positive A^, as in the Appendix, 
by the expression 

oo 

4iP,q,k) ^ d*{p,q,k):= Cn{p)Cn{q)Cn{k) (90) 

n=N+l 

oo 
n=N+l 

_|_ g-«2V"(p+cj"fc) _|_ gi2v^(-p+cj-fc) 

By the theorem A.l of the Appendix, the infinite oscillatory series converges 



as 



oo «QQ 

J2 n-^^/^e^^'^v^ ~ / y-'/^e'^'^^dy. (91) 

— at:, JN 



n=N+l 

The improper integral in (19T]) converges (except for A; 7^ 0) to 



X 



^-3/4g.2fcv^^^ ~ -{v^5(A;) + (l/4)|A;|-i/2 

1 -i)r (^-l/2,i2ViV|A;|) + {l + i)V i^-l/2,-i2^\k\^ 
-N-^/%\-^ s\Yi{2VN\k\) 

+i2-^\k\-^/\gn{k)[{l + i)T{l/2,i2VN\k\) (92) 
+ {l-i)T{l/2,-i2y/N\k\)]. 



Here, r(a,i2;) is the incomplete gamma function. Using the asymptotic 
representation of the incomplete gamma function for very large \/N\k\, 
we derive from ([90D dH]) ([92D that 

/■oo 

/ y-^/^e'^''^ dy = N'^^"^ cos{2VN\k\)5{k) + isgn{k) exp[i2VNk] + 0{l/N^/^\k\^) 

^ {N~'^/^\k\-^ [N-^/'^6{k) + isgn{k)] . (93) 
Therefore, by the equations we conclude that 

{2^fd*{j),q,k) ~ N-^'^{\p + q + k\-^[N-^l^5{p + q + k) + isgn{p + q + k)\ 

+ \q + k — \N~^/'^5{q + k — p) + isgn{q + k — p)~^ 
+ \q — k + p\^^ [N^^^'^6{q — k + p) + isgn{q — k + p)~\ 
+ \q — k — p\^^ [N^^^'^6{q — k — p) + isgn{q — k — p)~^} 
= ■.aNip,q,k). (94) 
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Thus, by equations (pO!) and (!93|) we obtain the asymptotic approximation 
(87r=^/2)353#,(p,q,k) (87r3/2)34^(p,q,k) 

3 



6=1 



hm 

{Ni,N2,Ns)^{oo,oo,oo) 



(96) 



6=1 
3 

= llW^'lPb + qb + h\-'[N;'^'6{pb, h) + 

6=1 

isgn(pfe + g;, + fcfe)] + .. + .. + ..}. (95) 

For completion of the step iii) in the strategy, we need to consider j products 
of the triple series as in (6.12). These multiple series obviously converge 
asymptotically to the products of terms as in [951 

As for the step iv) of the strategy, we notice that the analogous necessary 
condition 

3 

Yl^N,{Pb,qb,h) 

-6=1 

is satisfied by (p5!) . 

Now we have to investigate the most important step v) of the strategy. 
Suppose that for the S-matrix element (/ \ Sj \ i) in (1771) yielding the equa- 
tion (6.10), the corresponding Feynman graph involves the following specifi- 
cations: 

Ep =The number of external electron-positron or fermion lines, 
Eb =The number of external photon or boson lines, 
Ji? =The number of internal fermion lines, 
Ib =The number of internal boson lines, 
j =The number of corner or vertices. 

For a trilinear interaction as in the equation flSTl) . the following equations^ 
among various numbers hold: 

If = j-{Ep/2)>0, 
Ib = (j-^b)/2>0, 
J = If + {Ef/2)=2Ib + Eb>0. (97) 

The typical S-matrix element in (6.10), with the help of the modified Feyn- 
man rules in momentum space (Table-I), yields the asymptotic approxima- 
tion: 

<f\SjJi> ~ {const) / d^k(^)..d'^k(^ig) I d'^p(i)..d'^p(ip) 
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X {Y [7>{i)m + ^ ■■ [YP(If)u + ml] ^ 7p 

xo^Sv (P(i). -P(2), k) 6 - p(2)4 + .. 

c^Sv (-P(i-i), P0> k(j)) 5 (-po-i)4,P(i)4, %)4) I • (98) 



Integrating over all p^i)^ = —E(i) and k{i)^ = —^(i) variables we get from (l98i) 

<f\Sjj^\i> = {const.)6 (y2{E{i) + '^{i))) / rf^k(i)..c/^k(/^) / c/^p(i)..rf^p(/^) 

-1 



dk(j+i)^...dp(j+i)^... X <^ [7>(i)a + ••] •• [l^P{ip)b + 



[kli)k{i)c] ^ •• KB)%s)rf] ^ ^ ^'^\PWi ~Pm + 



Ni ^^^S {p(i)i - p(2)i + k(^i)i) + i sgn - p(2)« + A;(i)z) ) 
+similar terms . 

The Dirac delta function 5 ((^ {^{i) + '^(o)) indicates the conservation of 
total energy. The number of integration variables in (p9|) is 3 (Is + Ip) + 
(/^ + /f ~ J + 1) • At high energy-momentum, the factor with products of 
all |p|~^sgn(p) and |fc|~^sgn(fc) in the integrand, behaves asymp- 

totically as the power —Ip — 21 b — 3j. Therefore, counting the degree, k, 
of ultra-violet divergence with p''^,k~'^, sgn(A;) and |p|~^sgn(p) factors 
alone, is 



3{Ib + If) + {Ib + If-J + 1)-If-2Ib-3j = 1-Eb~{3/2)Ef. (100) 



The degree, k of divergence with product of all |p| ^S{p) and \k\ ^6{k) 
functions alone (after integrating out all (5-functions) is 

k := 4(Jb + If-j + 1)-If- 2Ib - 3j = 4 - 3j - Eb - (3/2)^^. (101) 

From (I loop and (1 10 II) we derive that 

K-K = 3(j - 1) > 0. (102) 

There exist other degrees of divergences arising from the mixed products of 
sgn(..) and These degrees are all less or equal to k. Therefore, the 

final criterion for "convergence" of the S-matrix element < f\Sjj^\i > or 
< f\Sj\i > is that 



K < 0, 
{3/2)Ef + Eb > 1. 



(103) 
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The above inequality proves that in case the number of external boson lines is 
at least two, the S-matrix elements "converge" in every order of j . However, 
in case the number of external bosonic lines is one, the S-matrix element 
vanishes by Furry's theorem. Therefore if there is at least one bosonic ex- 
ternal line, the S-matrix elements converge in every order. Now, the number 
of incoming external fermionic lines always equals to the number of outgoing 
fermionic lines. Therefore, we conclude from the inequality fll02p that in case 
the number of incoming external fermion lines is at least one, the S-matrix 
elements "converge" in every order j . Thus, our final conclusion is that the 
S-matrix elements in the discrete phase space formulation converges in every 
order provided there is at least one external line in the process. 

Let us consider some specific examples. In case of the self-energy of an 
electron, as discussed in section-5, the relevant numbers follows: 

Ef = 2, Eb = 0, k = -2< 0, 

{3/2)Ef + Eb = 3> I. (104) 

Therefore, by the equations fll02p . the S-matrix elements representing the 
self-energy of the electron "converges" in every order. 

As a second example, consider the self-energy of the photon. By the 
treatment in section-5, we deal with the numbers: 

Ef = 0, Eb = 2, k=-1<0, 

(3/2)Ef + Eb = 2>1. (105) 

Again, the corresponding S-matrix elements for the self-energy of the photon 
"converges" in every order. 

It can be shown that each of the so called primitive diagrams, "converges" 
in every order within the arena of the discrete phase space formulation of the 
S-matrix. 
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A Appendix: The distribution function (5*(p, q, k) 

We define the distribution function 5* by the following equation; 

oo 

6*{p,q,k) := ^en(p)a(g)a(fc), 

n=0 

(p,g,A;)GR^ (106) 



We assume that the right-hand side is at least summable in some sense. The 
distribution function 5* so defined, is totally symmetric with respect to the 
real variables p, q, k and thus satisfies 

6*{p, q, k) = 6*{q,p, k) = 6*{p, k, q) = 6*{k, q,p). (107) 

It has also the following properties by the consequences (A. 7) of paper-II: 

6*{p,q,k) = S* {p, -q, -k) = 6*{-p, -q, k) = 6*{-p, q, -k), (108) 
5#i-p, -g, -k) = 5*{p, q, -k) = 5*{-p, q, k) = 5*{p, -q, k), (109) 

Now we shall state and prove in a simple manner the following the- 
orem. 

Theorem A.l: The distribution function 5* is distinct from the Dirac 
(27r)(5 -distribution. 

Proof: Assuming the interchangeability of the summation and integration, 
and using the equations (II-A.I.l) and (II-A.I.19) we obtain 



« oo p « 

/ 'lM^*{p,q,k)dk = V / 



im{k)^n{k) dk 



in{p)in{q) 



= U{p)U{q)- (110) 

(Here the subscript m is not summed!) However, by the equation (II-A.I.16), 
we have 



R 



im{k)5{p + q + k)dk = ^m{-p - q) = ^mip + q) 

= (vr)^2-/2exp[(p-g)V2] 



m 

xE^^— ^■(^p)^^(^^)- (111) 

j=0 



Comparing the equations fill 01) and fill II) we conclude that S'^i^Pjqjk) ^ 
(27r)5(p + g + fc)n. 

Now we shall introduce another sequence of functions 



1/2^-1/4 



vr ' n 



Uk) : =Uk)ioinE{0,l}, 
cos{2k^/n) for n E {2,4,6,...}, 



i^T-^/^n-^/'^sm{2k^/E) for n G {3,5,7,...}. (112) 
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If we introduce a new inner-product and norm for a non-separable space® : / 
by 



# 



lim 



then it can be shown by flll2p that 

||e(k)-c(k)||^o 



2N+1 



n=0 



(113) 



(114) 



Therefore, it is reasonable to introduce an asymptotic approximation of 
(5*(p, q, k) in (11061) by the expression 



d*ip,q,k) : =J2Up)Uq)Uk) 

n=0 

00 

+ (l/87r3/2)^n-3/4 ^ 

n=2 

_|_ Q-i'^Vniq-k+p) _|_ ^i2^(q-k-p) 



^i2^{p+q+k) _|_ g-i2v^(ij+fc-p) 



;ii5) 



The convergence or summability of the above series is not obvious. However, 
we shall discuss some related results. We can cite^ the following theorem. 

Theorem A. 2: Let f and p be two real-valued and continuous func- 
tions defined over [a, 00). Let {&m}o° a monotone increasing sequence 
such that bo = a and limm^oo ~^ 00 • Moreover, let f{x) > and 
{—l)"^p{x) > for X e \brmbm+i\ C [a, 00) Then, the improper (Riemann or 
Lebesgue) integral f{x)p{x) dx converges to a iff the alternating series 



Eoo 
n= 







reference 9.) Now we shall prove the 



converges to a. (For proof see the 
bllowing corollary. 



Corollary A.l: There exists a sequence {Qnik)}^^ with the property 

[{n+e„(k)y, 

n+en{k)]l^ 



Qn{k) E (0, 1) for every n E N such that the series Yl'^=o 



converges to ^\a;')'i-^;3 Jq^ all k > and all [3 G (0, 1). 

Proof: It is known that for A; > and < /9 < 1, the improper integral 



sin(A;x) 



dx 



0+ 



X'' 



r(l-/j) sin[(l-/?)7r/2] 

{ky-f^ 

TT 

2(A;)i-/3cos[(l-/?)7r/2]r(/5)' 



(116) 
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If we choose a = 0, /(x) = x ^ , p{x) = sm{kx) {k > 0) , bn = {nn/k), and 
apply the preceding theorem-A.2, then we obtain that 



ni-/3) sin[(l-/3)7r/2] _ >^ 

(^)l-/3 2^ 



sinffcx) , 
— ax 



;ii7) 



n=0 

Applying the mean- value theorem for an integral, we conclude that 

y= [''+e„wp to.=o,nee„w.(o,i).n (iis) 

We can generalize the preceding corollary for any k ^ and (3 G 
(0, 1) to the following equations: 

^ sin[(n + 9n(A;))7r] _ sin(fcx) 



sgn(A;)r(l-/3) sin[(l - /3)7r/2] ^ ^^^^^ 



|A;|i-/3 



cos[(n + 6„(/c))7r] cos{kx) ^ 



cos[(l - /3)V2]{^]^T^ - 2r(-/3)|A;|'^5(fc)}, 



y. exp{z[(n + e#(fc))7r]} ^ H dx 

= ^|^^exp{^sgn(fc)(l-/3)7r/2}(120) 

-2cos[(l -/3)7r/2]r(-/3)|A;|^5(A;). 

Here, is the Dirac delta distribution. 

It is very plausible from flllQl ii) that the series d*{p,q,k) in (11151) 
converge. In fact, from (11191 . ii, iii) we shall choose an approximation for 
q, k) by the following relation: 

8d*{p, g, k) ^ vr-^/^ / y'^^^ (exp{i[2(p + q + k)^} 

+ exp{-2[2(g + k-p)y/y]} + exp{-i[2(g - k+p)y/y]} + exp{i[2(g - k-p)y/y]}) dy 

POO 

= V2(7r)-3/2 / x-^/^{exp[i{p + q + k)x] 
Jo+ 

+ exp[— ^(g + k — p)x] + exp[— z(g — k + p)x] + exp[z(g — k — p)x]}dx. (121) 
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Now, choosing the special case /? = 1/2 and recalhng r(l/2) = 
v/vr, r(-l/2) = -2^/?, we obtain from (A.12iii) 



r x-'/\^^^dx = V^{2^\6{k) + il±ll^}. (122) 
Jo+ \/2|A;| 



Substituting (11221) into (11211) . we finally obtain an asymptotically "closed 
form" : 



2Tid*{p,q,k) = {a/Ip + q + k\S{p + q + k) + ^/\q + k- p\S{q + k - p) 
+ V\q - k + p\S{q - k + p) + ^/\q - k - p\S{q - k - p)}. 
[1 + i sgn(p + q + k)] [1 — « sgn(g + k — p)] 



+ imv 



y^\p + q + k\ \^\q + k-p\ 



^ sgn{q-k + p)] ^ [1 + i sgn(g - A; - p)] ^ ^^^^^ 

^y\q -k + p\ -k-p\ 

We shall now state and sketch briefiy the proof of the remarkable 
properties of the distribution function d'^ip^q^k) in (I123p . 

Theorem A. 3: The distribution function d'^{p,q,k) satisfies exactly the 
symmetry properties (A I1071 [T08| I109p of d'^ip^qjk). 

Proof: The proof is straight forward (using I123p and is skipped. 
We define the three dimensional distribution function 

3 

4(p,q,k) ■.= \{d*{p,,q,,k,). (124) 

i=i 

In section- V, we shall use the asymptotic approximation of the 5f- 
distribution function by putting 

53#(p,q,k)^4(p,q,k). (125) 
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